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ABSTRACT 

We use high-resolution TV-body numerical simulations to study the number of pre- 
dicted large-separation multiply- imaged systems produced by clusters of galaxies in 
the SDSS photometric and spectroscopic quasar samples. We incorporate the conden- 
sation of baryons at the centre of clusters by (artificially) adding a brightest central 
galaxy (BCG) as a truncated isothermal sphere. We make predictions in two flat cos- 
mological models: a ACDM model with a matter density £! m .o = 0.3, and <r$ = 0.9 
(ACDM0), and a model favoured by the WMAP three-year data with f2 m;0 = 0.238, 
and cxs = 0.74 (WMAP3). We found that the predicted multiply-imaged quasars with 
separation > 10" is about 6.2 and 2.6 for the SDSS photometric (with an effective 
area 8000 deg 2 ) and spectroscopic (with an effective area 5000 deg 2 ) quasar samples 
respectively in the ACDM0 model; the predicted numbers of large-separation lensed 
quasars agree well with the observations. These numbers are reduced by a factor of 
7 or more in the WMAP3 model, and are consistent with data at < 8% level. The 
predicted cluster lens redshift peaks around redshift 0.5, and 90% are between 0.2 and 
1. The ratio of systems with at least four image systems (JVi, ng > 4) and those with 
N img > 2 is about 1/3.5 for both the ACDM0 and WMAP3 models, and for both the 
photometric and spectroscopic quasar samples. We find that the BCG creates a central 
circular region, comparable to the Einstein ring of the BCG, where the central image 
disappears in the usual three-image and five- image configurations. If we include four 
image systems as an extreme case of five- image systems (with an infinitely demagnified 
central image), we find that 68% of the central images are fainter by a factor of 100 
than the brightest image, and about 80% are within 1.5" of the BCG. 

Key words: gravitational lensing - galaxies: clusters: general - cosmological param- 
eters - dark matter 



1 INTRODUCTION 

The number of multiply-imaged quasars lensed by galaxies 
has now reached roughly one hundred^. The typical separa- 
tion of these systems ranges from 0.3 — 6". They provide an 
valuable sample to constrain the cosmological constant, the 
Hubble constant, and the mass profiles of lensing galaxies at 
intermediate redshift (for an extensive review, see Kochanek, 
Schneider & Wambsganss 2006 and references therein). 

The search for multiply-imaged quasars by clusters of 
galaxies has been less successful. The initial search for wide- 
separation radio sources yielded no successful candidates be- 



* E-mail: lglSshao .ac.cn 

1 http:/ /cfa-www. harvard.edu/castles/ 



tween 6 - 15" JPhillips et al.l l2001al). a nd from 15 - 60" 
(|Phillips et al.ll2001bh . lOfek et all (|200l"h also failed to find 
any wide-separation quasars between 5 — 30" in the FIRST 
20-cm radio survey. The lack of cluster lensed radio sources is 
likely due to the small number of radio sources surveyed (~ 
20,000). The breakthrough came from the Sloan Digital Sky 
Survey (SDSS) where a large number of optical quasars be- 
came available. So far, there are approximately 46,420 spec- 
troscopically confirmed broad-line quasars in SDSS Date Re- 
lease 3. Two spectacular cluster lenses were discovered in 
the SDSS (not limited to DR3, see H4.2p . The first case was 
SDSS J1004+4112, a quadruply-imaged system with an sep - 
aration of 14.6 arcsec (jlnada et al.ll2003l ; fOguri et al.ll2004r ). 
The cluster is at redshift of 0.68, and the lensed quasar is 
at redshift 1.734. A faint fifth image was later discovered 
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l|lnada et all 120051 ). as were many othe r multiply-imaged 
background galaxies (jSharon et al.ll2005M . The second case, 
SDSS J1029+2623, was a double-image system with a sepa- 
ratio n of 22.5 arcsec. T he lensing cluster is likely at redshift 
0.55 (|lnadaet al.ll2006T ). and the background quasar at red- 
shift z s = 2.197. 

Cluster lenses, just like galaxy-scale lenses, provide 
valuable constraints on the lens mass profiles. They are also 
important because the number of such lenses depend quite 
sensitively on the cosmology, in particular the matter power- 
spectrum normalisation, as, and the matter density, f2 m ,o. 
This sensitivity has already been explored using giant arc 
statistics where the lensed sou rces are background gala x- 
ies, rather than 
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quasars (e.g., LM cncghct ti et al.l l2 003a 
L l2004l:lTorri et al .1120041: iDafal et al 11200 



Horesh et alj 12003: iLi et all 120051 ; lOguri et al l \2003 ). The 
analysis by Li et al.l ( 2006bl ) prefers a high as universe. How- 
ever, the study and other similar ones suffer from one de- 
ficiency, namely the properties of the lensed background 
galaxy population, in particular their redshift distribution, 
is not well known, so the conclusion is somewhat uncertain. 

In comparison, the quasar population has a much 
better-known redshift distribution, and so does not suffer 
as much from the same deficiency, hence a careful analysis 
of the statistics of multiply-imaged quasars is complemen- 
tary to those of giant arcs. Of course, the number of known 
multiply-imaged quasars so far is very s mall (two) , compared 
to hundreds of known giant arcs (e.g.. iLuppino et ail 19991: 
Gladders et~aill2003l: IZaritskv fe Gonzalezll2003l : ISand et ail 



20051 : iHennawi et alj|2006l ). Nevertheless, even this small 
sample pro vides interest i ng co nstraints. For example, the 
analysis of lOguri et alj l|2004l ) of the first cluster lens, 
J1004+4112, concluded that the as value must be high, 
as — lto'2 (95% confidence) even for an inner density pro- 
file, p cc r~ 1,5 , considerably steeper than the NFW pro- 
file (p oc r _1 for small r, Navarro, Frenk & White 1997). 
Many stu dies of clusters of galaxies prefer the NFW pro- 



file (e.g., van de Marel et "aiTl2000l : IComerford et~ail 120061 : 



IVoigt fc Fabianll2006l ). For such a shallower profile, the re- 
quired as would be even higher. 

Another advantage of lensed quasars concerns the lens 
selection. For giant arc surveys, observers first select clusters 
from optical or X-ray data and then take deep followup im- 
ages, hence the discovery of giant arcs depends on how clus- 
ters are selected in the first place. On the other hand, lensed 
quasars are usually identified by examining quasars to see 
whether they are multiply imaged. In other words, current 
arc surveys only observe biased directions whereas quasar 
lenses survey random line-of-sight. In addition, quasars are 
point sources, so one does not have complications arising 
from seeings, background source size and ellipticity distri- 
butions etc, thus it is much simpler to model the quasar 
source population. 

There have been many previous works to predict the 
abundance of large-separation quasar l ensed by massive 
dark halos using simple spherical models l|Maoz et al 19971: 
Sarbu et ail l200ll: iKeeton fc Madaul " 
Li & Ostriker 2002; Rusin 2002; 
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200 ll; IWvithe et al l 
ILi fc Ostrikerl | 2003 | : 
Qguril 120031: iHuterer fc Mai |2004|; iLopes fc MiiieTbooi 
Kuhle n et al. boO^IOguri" et al . 2004; Chen 200l)- In partii 



Table 1. Parameters for the two cosmological numerical simula- 
tions. The columns are cosmology, box size, number and mass of 
dark matter particles, and the softening length. 



model 


box size 


AW 


rriDM 


softening 




(h- 1 Mpc) 






(h- 1 kpc) 


ACDM0 


300 


512 3 


1.67 x 10 10 


30 


WMAP3 


300 


512 3 


1.32 x 10 10 


30 



ular, ILopes fc Miner] l|2004h . using spherical clusters, found 
that wide-separation lensed quasars can be a powerful 



tool to measure ag an d £l m ,o- Using triaxial halo models, 
lOguri fc Keetonl l|2004h computed the abundance of quasars 
multiply imaged by clusters analytically and found that 
triaxia lity significantly incre ases the expected number of 
lenses. [Hennawi et alj (|2007bl) presented the first predictions 
of multiply-imaged quasars using N-body simulated clusters 
in the a flat cosmology with O m ,o = 0.3, Qa,o = 0.7 with 
as = 0.95. In this paper, we compare the predictions for the 
cosmology favoured by the three-year WMAP result with 
the usual concordance cosmology. We also make more de- 
tailed predictions for cluster lenses as a function of image 
multiplicity, and study the brightness of the central images 
and their offset from cluster centres. 

The paper is structured as follows. In §2, we discuss 
the simulations we use and our lensing methodology. In §3, 
we discuss the background source population, namely the 
quasar luminosity function. In §4, we present the main re- 
sults of the paper, and we finish in §5 with a summary and 
a discussion of the implications of our results. 



2 NUMERICAL SIMULATIONS AND 
LENSING METHODOLOGY 

In this paper, we utilise two sets of simulations, one for the 
'standard' ACDM model with as = 0.9, and the other for the 
cosmological model given by the recent three-year WMAP 
data with a 8 = 0.74 (|Spergel et al.ll2006l) . For brevity, these 
two models will be referred to as ACDM0, and WMAP3, 
respectively, and the corresponding cosmological parameters 
are: 

(i) ACDM0: n m , = 0.3,fi A ,o = 0.7, h = 0.7, a 8 = 
0.9, n = 1; 

(ii) WMAP3: fi m , = 0.238, « A ,o = 0.762, h = 0.73, a 8 = 
0.74, n = 0.95, 

where h is the Hubble constant in units of 
100 km s" 1 Mpc" 1 , and n is the spectral index of the power 
spectrum. The assumed initial trans fer function in each 
mode l was generated with CMBFAST (|Seliak fc Zaldarriagal 
1996). Notice that both as and f2 m ,o differ in these two mod- 
els, and both will impact on the number of multiply-imaged 
quasars. 

We us e a, vectorised-parallel P 3 M code 

ilJing fc Sutol 120021) and a PM-TREE code - GAD - 
GET2 (|Springel. Yoshida. fc White I l200ll ; ISpringel 1 120051 ) 
to simulate the structure formation in these models; the 
details of the simulations are given in Table 1. Both are 
A^-body simulations which evolve JVdm = 512 3 dark matter 
particles in a large cubic box with sidelength equal to 
300/t" 1 Mpc. The ACDM0 simulation was performed using 
the P 3 M code of Jing & Suto (2002) and has been used 
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m iLi et all (|2005l ) to study the properties of giant arcs. 
We refer the readers to that paper for more details. The 
WMAP3 simulat ion was carried out with GADGET2, and 
has been used bv lLi et alJ (|2006bh to compare the statistics 
of giant arcs in the ACDMO and WMAP3 models. Similarly, 
in this work we will use them to compare the statistics of 
multiply-imaged quasars in these two cosmological models. 
We return to the issue of cosmic variance in i|4.4l 

Dark matter halos are identified with the friends-of- 
friends method using a linking length equal to 0.2 times 
the mean particle separation. The halo mass M is defined 
as the virial mass enclosed within the virial radius accord- 
ing to the spherical collapse model jKitavama fc Sutoll 19961 ; 
iBrvan fc Normanlll998l ; Ijing fc Sutoll2002l) . 

For any given cluster, we calcu late the smoothe d surface 
density maps using the method of ILi et all l|2006al '). Specif- 
ically, for any line of sight, we obtain the surface density 
on a 1024 x 1024 grid covering a square of (comoving) side 
length of 4/i -1 Mpc centreed on each cluster. The projection 
depth is also chosen to be 4fo _1 Mpc. Notice that the size 
of the region is large enough to include all particles within 
the virial radius (in fact, several virial radii for small clus- 
ters). Particles outside this cube and large-scale structures 
do n ot contribute significantly to the le nsing cross-section 
(e.g. ILi et alj|200a : iHennawi et alj|2007al ). so we will ignore 
their contributions here. 

Our projection and smoothing method uses a smoothed 
particle hydrodynamics (SPH) kernel to distribute the par- 
ticle mass on a three-dimension grid and then the surface 
density is obta ined by integrating along the line of sight (see 
ILi et al.ll2006al for more detail) . In this work, the number of 
neighbours used in the SPH smoothing kernel is fixed to be 
32. We calculate the surface density along three orthogonal 
directions for each cluster. Once we obtain the surface den- 
sity of a dark matter halo, its lensing potential, 4>dm, can 
then be easily calculated using the Fast Fourier Transform 
(FFT). 

Many clusters of galaxies have a very luminous central 
galaxy. The effect of s uch brightest cluster ga laxies (BCGs) 
has been quantified bv lHennawi et al.l(|2007bl ). For multiply- 
imaged systems with separations < 20", BCGs can en- 
hance t he cross-section by ~ 50%. Following IHennawi et al.l 
(2007b), we also artficially add a BCG at the centre of each 
cluster. This is obviously a simplification as BCGs are as- 
sembled as a function of time; we briefly return to this is- 
sue in the discussion. The BCG is modelled as a truncated 
isothermal sphere. The mass of BCG is taken to be a fixed 
fraction of the host halo, Mbcg = 3 x 10 _3 Mfof and its 
velocity dispersion is given by 

2/15 

(1) 



300 kms- 



Mfof 



10 15 M 



Th e truncation radius is g iven by r max = GA/bcg/(2<7 2 ). 
SeeE cnnawi et alJ (|2007bh for more details. 

In the following we will only present results includ- 
ing such a BCG, which approximately accounts for the 
effects of baryonic cooling at the centre of clusters. The 
effect of baryons in non-central galaxies is likely to be 
smaller. As the BCG is modelled as a truncated isother- 
mal sphere, its lensing potential, 0bcg, can be calculated 
analytically. The total lensing potential is a sum of the con- 
tributions from the (pure) dark matter halo and the BCG, 



4> — 0dm + 0bcg. The deflection angle (a) and magnification 
(fx) can be obtained by taking first-order and second-order 
derivatives with respect to x and y, a x = (j> x , a y = <j)y, and 
(j, = [(1 - <f> xx )(l - <t>yy) ~ (Ply]' 1 - We calculate the BCG 
contribution to the deflection angle and potential analyti- 
cally, and the dark matter halo contribution using FFT as 
discussed above. 

Using these derivatives, we can create the mapping from 
the image plane to the source plane and obtain the critical 
curves and caustics. Critical curves are a set of image posi- 
tions with jj, = oo, and the caustics are the corresponding 
source positions. Caustics play a crucial role in gravitational 
lensing as they divide the source plane into regions of differ- 
ent image multiplicities. Whenever a source crosses a caus- 
tic, the image number either increases or decreases by two. 
For a quasar to be multiply-imaged, it must lie inside one of 
the caustics. 

To calculate the cross-section of a given separation of a 
multiply imaged quasar, we first find the smallest rectangle 
which encloses all the caustics, and then a regular grid is cre- 
ated dividing this rectangle into pixels of 0.3" by 0.3". For 
each source position located on the grid, we find the corre- 
sponding image positions and magnifications (p) using the 
Newton-Raphson method. For a given multiple-image sys- 
tem, the number of images one can observe, Ni ms , depends 
on the magnification of each image, the source intrinsic mag- 
nitude and the survey magnitude limit. For such systems, we 
sort the images according to the image magnifications and 
record the number of images -/Vi mg and the maximum sep- 
aration A9 between these iVi mg images. We then examine 
whether the multiple images satisfy the observational selec- 
tion criteria (see jj3.1|) . For given lens and source redshifts, 
z\ and z a , collecting systems that pass the selection allows us 
to construct a cross-section cr(> M mg , > 9\z\, z s ) for lensing 
systems with at least iVimg images and separations larger 
than 6. 

The SDSS photometric and spectroscopic quasar sam- 
ples are magnitude limited, multiply-imaged quasars will be 
over-re presented in s uch samples due to the lensing magnifi- 
cation (Turner 1980). To account for this magnification (or 
amplification) bias, we further divide the cross-section into 
50 bins of Am = — 2.5 log from —5 mag to 5 mag with a 
bin size of 0.2 mag, and obtain a differential cross-section 
da/dAm(> A r i mg , > 9\Am, zi, z s ). 

As in ILi et al. I l|2005h . we can obtain the total cross- 
section per unit comoving volume by summing over the con- 
tributions of all the clusters: 



&(> Aim g > > 6\Am, zi, z s ) = 



cr t (> Aimg, > 9\Am, z u 
V 



',(2) 



where <Tj(> JVi mg , > 9\Am, z\, z s ) is the average cross-section 
of the three projections of the i-th cluster at redshift zi, z s 
is the source redshift, and V is the comoving volume of the 
simulation box. Then, accounting for the magnification bias 
(for more detailed descriptions, see Sect. 3.1), the predicted 
number of large-separation lens can be derived as: 



N(>N img ,>9) 



Jo" dz - fo S dz i /*» dAm 
^(>N img ,>9\Am,z h z s ) 



x n(< i — Am, z s ) (1 + z\) 



3 dVp(fl) 



(3) 
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where n(< i, z s ) X dz s is the number of quasars brighter than 
the SDSS i-band magnitude i in the redshift interval from 
z s to z s + dz s , io is the flux limit, and dV p (z\) is the proper 
volume of a spherical shell with redshift from z\ to z\ + dz\. 
We used 22 simulation output with redshift from 0.1 to 2, 
and selected the 600 (for ACDMO) or 400 (for WMAP3) 
most massive clusters in each output to calculate the total 
lensing cross-sections. The lowest masses of clusters in these 
two cosmologies are 7.8 x 10 13 /i _1 A/ Q and 5 x 10 13 /i _1 M© 
respectively at redshift 0.5. The integration step size of z\ is 
the same as the redshift interval of simulation output {5z\ ~ 
0.1) and the source redshift is at 0.7, 1.0, 1.5, 2.0, 2.5, 3.0, 
3.5, 4.0, 4.5, 5.0 and 5.5. 



3 THE QUASAR LUMINOSITY FUNCTION 

To make accurate predictions for the frequency of multiply- 
imaged quasars, we must specify the quasar luminosity func- 
tion as a functio n of redshift and luminosity. Following 
iBovle et al.l (|2000h . we will use a double power-law to pa- 
rameterise the quasar luminosity function 

$*(z) 

<f>(M,z) = 1()0 . 4(/ 3 i+1)[M _ M . {z)] - 10 o.4(^ h +l)[A/-Af«(z)] ' ( 4 ) 

where M is the absolute magnitude in a filter, and $* and 
M* axe the characteristic density and magnitudes respec- 
tively in a certain band. 

For the low redshift (z < 2.1) quasars, we use the pa- 
rame ters from the 2dF-Sp SS LRG and QSO (2SLAQ) Sur- 
vey |Richards et alj 120051 ) . The luminosity function is fit- 
ted in the SDSS fir-band with fa = -1.45, (i h = -3.31 and 
$* = 1.83 x 10~ 6 Mpc~ 3 mag _1 . The characteristic lumi- 
nosity in the g-band, Mg, is described by a second-order 
polynomial of redshift: 

M* g (z) = M 9 *(0) - 2.5(fciz + k 2 z 2 ), (5) 

where M*(0) = -21.61, ki = 1.39 and k 2 = -0.29. 

For high redshift, z > 3.5, quasars, we take the parame- 
ters from lFontanot et all (|2006T l. which were obtained from a 
joint analysis of the GOODS and SDSS data. The best fit to 
the luminosity function as a function of the monochromatic 
1450 A magnitude, Mas, gives fa = -1.71, p h = -3.31; the 
data also prefers a pure density evolution with an exponen- 
tial form: 

$*(z) = $(_ 2)e fc *'< 1+z >- 3 l, (6) 

where k z = -1.37 and $* z=2 ) = 1-67 x 10" 6 Mpc" 3 mag -1 . 
The characteristic monochromatic luminosity, M* 45 , is fixed 
at M 145 = -26.43. 

For a given flux limit io, we can easily obtain the total 
number density of quasars n(< io,z) in the redshift space 
for z < 2.1 and z > 3.5. However, the quasar luminosity 
function is not well constrained from redshift 2.1 to red- 
shift 3.5. For a given magnitude limit io, we estimate the 
number density of quasar with redshift 2.1 < z < 3.5 by 
linear interpolation as a function of redshift of the quasar 
number densities with the same absolute magnitude. For a 
given redshift z, we first calculate the corresponding abso- 
lute magnitude, Mf(z) for each apparent magnitude i. We 
then convert this, using k-corrections (see i)3.1fl . into corre- 
sponding g-band and the 1450 A monochromatic absolute 
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Figure 1. The left panel shows the redshift distribution of 
quasars for a flux limit of i < 22.7. Three redshift intervals 
(z s ~ 1,2 and 4) are indicated by the vertical lines. The shaded 
region corresponds to the 60% completeness we assumed. The 
right panel shows the surface density of quasars as a function of 
the SDSS i-band magnitude for the three source redshift intervals. 
The solid and dashed lines show the shapes used by us and by 
iHennawi et al. (2007b) respectively. They are in good agreement 
for the two low redshift intervals, but for high reshift interval, our 
amplitude is higher by roughly a factor of 2 than theirs between 
i = 18 and 22. 

luminosities at Mg and Mi 45 at redshift z. Since we know 
3>(M°,z = 2.1), $(M? 45 ,z = 3.5), we obtain the $(M;,z) 
by linear interpolation: 

_ $(M°,2.1) x (3.5-z)+$(Mi°45,3.5) x (z - 2.1) 
1 t,Z> 3.5-2.1 

We can then simply integrate the above equation to derive 
the number density of quasars per square degree for any 
limiting magnitude i. 

We also consider the incompletenes s of qua sars at red- 
shift from 2.3 to 3.6 as in lHennawi et al.1 l |2007bl ). The SDSS 
spectroscopic and photometric surveys depend on colour se- 
lection and the colours of quasars and stars are not easy to 
disentangle in this redshift range. So we assume that only 
60% of quasars in this redshift interval are found by the 
SDSS. 

3.1 The cross filter k-correction, quasar samples 
and lens selection 

The SDSS spectroscopic and photometric survey have flux 
limits in the i-band but the parameters in the above lumi- 
nosity functions (in eqs.HKJ) are given in different bands, M g 
and M145 . To apply these quasar luminosity func tions to the 
i-ban d, we have to do cross filter k-corrections dHogg et al.l 
2002). We use the results of I Richards et all (|2006h who pro- 
vided the k-correction from the SDSS i-band to the SDSS 
g-band and the monochromatic 1450 A absolute magnitude 
assuming a power-law qua s ar spec tral index, a v = —0.5. 

As in IHennawi et al.l ||2007bJ ), we use the same pho- 
tometric and spectroscopic quasar samples. For the SDSS 
spectroscopic quasar sample, the covered area is about 5000 
deg 2 and the magnitude limit is % < 19.1 for low redshift 
(z < 2.5 ) quasars and i < 20.2 for high redshift (z > 2.5) 
quasars ((Schneider et al. II2005I ). The faint SDSS photomet- 
ric quasar survey covers a larger area, ~ 8000 deg 2 , with a 
flux limit of i < 21.0 for z < 2.5 and i < 20.5 for z > 2.5 
l|Richards et al.ll2004D . To select lenses in the photometric 
and spectroscopic quasar samples, we require the brightest 
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image to be brighter than the magnitude limit given above, 
and other images to be brighter than i = 21.0, after account- 
ing for their magnifications. The fainter limit for other im- 
ages is because of fiber collisions in the SDSS, usually only 
one quasar image will be spectroscopically targeted, while 
other quasar candidates have to be confirmed with other, 
typically larger telescopes that can reach much fainter mag- 
nitudes. This is what actually occurred for the discovery of 
the first two lens es, and so we w i ll adop t this selection cri- 
terion, similar to lHennawi et all (|2007bl) . Notice that some 
images may still be missed if they are below the flux limit, 
even after the magnification, so the observed number of im- 
ages may be smaller than the intrinsic number. 

Notice that the photometric sample goes almost two 
magnitudes deeper than the spectroscopic sample for 
quasars below redshift 2.5 (magnitude limit 21.0 vs. 19.1), 
but only slightly deeper for high redshift quasars (20.2 vs. 
20.5). As a result, the photometric sample has many more 
low-redshift quasars compared with the spectroscopic sam- 
ple. This will be reflected in the source and lens redshift 
distributions for multiply-imaged quasars (see Figs. [S] and 

Fig. [1] shows the redshift distribution of quasars for 
i < 22.7. It shows the well known peak around redshift 2, 
and the sharp decline beyond redshift 4. The shaded region 
corresponds to the 60% completeness we assumed. This re- 
sult is in good agreem ent with the left panel of Fig. 1 in 
iHennawi et al.l (|2007bh . Notice that the geometries in the 
ACDM0 and WMAP3 cosmologies are different, due to the 
small difference in the matter density. However, even for a 
source at redshift 5, the angular diameter distance differs 
only by < 8%. Furthermore, we reproduce the quasars dis- 
tribution as a function of redshift and magnitude (see Fig.[T]) 
using our luminoisty functions. This distribution is directly 
comparable with the observations and model independent. 
Our calculation in each cosmology is based on this quasar 
distribution. So the small cosmological dependence of the 
luminosity function is accounted for. The left panel also in- 
dicates three redshift intervals {z ~ 1, 2 and 4), for which we 
show the surface density of quasars as a function of magni- 
tu de limit. The dashed lin es show the corresponding results 
ofl cnnawi et al l (|2007bh . Our results are consistent with 
theirs for the two low redshift intervals, but for the high- 
est redshift interval (2 ~ 4), our q uasar surface d e nsity is 
roughly a factor of two of that of IHennawi et all |2007b) 
between i = 18 and 22. As we use the updated results by 
Rich ards et al. (2005), our luminosity funct ion may be more 
realistic than that in IHennawi et~ai1 (|2007bh . It is worth em- 
phasising, however, that this difference has only moderate 
effects on the number of predicted multiply-imaged quasars 
because most (~ 75% and ~ 70%) of the lensed sources 
are below redshift 3 for the photometric and spectroscopic 
samples. 



4 RESULTS 

Armed with knowledge of the background source popula- 
tion and cluster lensing cross-sections, we can now make de- 
tailed predictions for multiply-imaged quasars in the SDSS. 
We start with an example cluster and illustrate possible 
multiply-imaged configurations in £14. 1 1 We then discuss sta- 
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Figure 2. The surface density map (in units of the critical surface 
density) for the most massive cluster at redshift zj = 0.5 in the 
ACDM0 simulation (in one projection). The cluster virial mass 
is 1 X 10 lr 'h~ 1 Mq. The colour bar at the right shows the scales. 
Notice that the main cluster (located at the origin) is merging 
with a sub-cluster at the left. The critical curves (appropriate for 
a source at redshift 2.5) are shown as the green curves (see also 
Fig. E]) 



tistical properties of multiply-imaged quasars, including the 
total number and partitions in different image multiplicities, 
for both the ACDM0 and WMAP3 cosmologies in £2] In 
i |4.3l we examine the properties of central images, and finish 
with a discussion of cosmic variance in <j4 .41 



4.1 



An example cluster: caustics, critical curves 
and images 



Fig. [5] shows the projected surface density map (along one 
direction) for the most massive cluster in the ACDM0 cos- 
mology at redshift 0.5. This cluster has a mass of about 
10 15 /i _1 M . It can be seen that in addition to the mas- 
sive cluster located close to the origin, there is another 
sub-cluster located roughly 40" away (corresponding to 
170/i _1 kpc away in projection) to the left merging with the 
main cluster shown at the origin. 

The critical curves and caustics for this cluster are 
shown in Fig. [3] Both are quite elongated along the hor- 
izontal axis, partly due to the shear of the merging sub- 
cluster. Notice that there is also a small, inner critical curve 
associated with the sub-cluster. The caustics consist of an 
inner ellipse and an outer diamond (there are small intersec- 
tions, however). There are also high-order singularities such 
as swallowtails in the caustics. Sources located inside such 
high-order singularity regions can produce image numbers 
higher than 5 (see the bottom right panel). 

The top left panel of Fig. [3] shows an image configura- 
tion where the source is located just inside the outer dia- 
mond caustic. All three formed images are quite far away 
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Figure 3. Example image configurations for the most massive 
cluster at Z\ = 0.5 in our ACDM0 simulation shown in Fig. [2] The 
source is assumed to be at redshift z s = 2.5. The cluster mass is 
1 X 10 15 /i~ 1 Mq. The blue and black curves are the caustics and 
critical curves respectively. For each panel, the red cross shows 
the source position and the corresponding images are shown as 
triangles. The bottom left panel shows an additional source po- 
sition (green cross) and the corresponding images are shown as 
green triangles. When a source lies in the shaded circle, the central 
image disappears. This is a result of approximating the Brightest 
Cluster Galaxy (BCG) as a truncated isothermal sphere (see the 
discussion in i|4,lj l. 



from the centre of the main cluster. When the source moves 
into the inner ellipse, we normally expect five images. How- 
ever, as a consequence of modelling the BCG as a truncated 
singular isothermal sphere, we find that there is a roughly 
circular region (shaded in the Figure) within which the cen- 
tral, fifth image disappears - the central image has been 
'swallowed' by the BCG. The radius of this region is roughly 
one angular Einstein radius corresponding to the isothermal 
sphere, 4-k(<t/c) 2 D\ b / D b , where D\ B is the angular diameter 
distance between the lens and source, and a is the velocity 
dispersion (cf. eq. [TJ. Such a four-image configuration is il- 
lustrated in the top right panel. We can regard such cases 
as an extreme example of a five-image configuration where 
the central image has been infinitely de-magnified. 

The bottom left panel in Fig. [3] shows a source position 
that produces the usual five-image systems. As expected, the 
fifth image is close to the cluster centre. There is, however, a 
different kind of five image configuration, which is produced 
by the presence of the sub-cluster. As we only put in (by 
hand) a BCG at the centre of the main cluster, the central 
image induced by such sub-clusters is always present, and 
will be located close to the centre of the sub-cluster, and in 
general much further away from the main cluster centre (see 
also Fig.©. 

For rare cases, when the source is located in regions of 
high-order singularity, the image number can exceed 5. One 
example is shown in the bottom right panel. For this case, 
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Figure 4. The top panel shows the number of predicted large- 
separation multiply-imaged quasars as function of separation for 
the ACDM0 (top two curves) and WMAP3 (bottom curves) mod- 
els, while the solid and dashed lines are for the SDSS photometric 
(with an effective area of ~ 8000 deg 2 ) and spectroscopic (with an 
effective area of ~ 5000 deg 2 ) quasars sample respectively. The 
bottom panel shows the ratio of the predicted numbers in the 
two cosmologies. Notice that the ratio increases as the separation 
increases. 



four images are clustered close together, tangentially aligned 
with a critical curve, while the other three images are fur- 
ther apart, ft can be seen from Fig. [3] that such high-order 
singularity regions are rare, and so we do not expect a sig- 
nificant fraction of sources to have image numbers exceeding 
five. 

For other clusters in our simulations, the image con- 
figurations are similar except that for some small clusters, 
we can also have two-image configurations. These are the 
usual three-image systems for which the central image has 
been 'swallowed' by the BCG, similar to the four-image case 
discussed above. 



4.2 Expected number of multiply-imaged quasars 
in the SDSS 

The expected number of multiply-imaged quasars in the 
SDSS is shown in Fig. 2] for the photometric (solid) and 
spectroscopic (dashed) samples, in both the ACDMO (top 
two curves) and WM AP3 (bottom curves) cosmologies. In 
the ACDMO cosmology, we expect 6.2 and 2.6 lenses with 
separations larger than fO" for the photometric (with an 
effective area of ~ 8000 deg 2 ) and spectroscopic (with an 
effective area of ~ 5000 deg 2 ) samples respectively. The cor- 
responding numbers are much lower in the WMAP3 model, 
about 0.85 and 0.36 respectively. The bottom panel of Fig. 
|3] shows the relative ratios in these two cosmologies - the 
ratio is about ~ 7 at fO", but increases to about ~ ff at 
30", for both samples. 

The large reduction factor in the WMAP3 model is a 
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Figure 5. The number of large-separation multiply-imaged 
quasars as function of separation and image multiplicity for the 
SDSS photometric (with an effective area of ~ 8000 deg 2 ) and 
spectroscopic (with an effective area of ~ 5000 deg 2 ) quasar sam- 
ples. The top and bottom panels show the results for the ACDMO 
and WMAP3 cosmologies respectively. For each panel, from top 
to bottom, the lines are for N- lIns > 2, 3, 4, 5 respectively. The 
two crosses show the predicted positions of the two known large- 
separation lenses. 
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Figure 6. The lens redshift distribution in the ACDMO cos- 
mology (left) and WMAP3 cosmology (right). The area under 
each curve is normalised to the total number of multiply-imaged 
quasars. The two vertical bars show the lens redshifts of SDSS 
J1004+4112 (0.68) and SDSS J1029+2623 (0.55). The shapes are 
similar in these two cosmologies. Notice that the median lens red- 
shift is slightly lower for the photometric quasar sample (see £[3j) - 
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direct result of the lower abundance of clusters in this cos- 
mology compared with ACDMO (see Fig. I10|l because of the 
lower as and lower mass density f2 m ,o. The reduction factor 
is larger for more massive clusters. This explains the separa- 
tion dependence of the ratio because larger separation lenses 
are produced by more massive clusters. The reduction factor 
is comparable to that of giant arcs (~ 6) for 9 > 10" but 
larger by a factor of 2 for 9 > 30". The substantially lower 
number of predicted multiply-images in the WMAP3 model 
may be somewhat difficult to reconcile with the data; we 
will return to this point at the end of this section. 

We can further divide the predicted number of multiply- 
imaged quasars by clusters according to the image multi- 
plicity. The results are shown in Fig. [5] for the two quasar 
samples in the ACDMO and WMAP3 cosmologies. While the 
amplitudes differ, the relative ratios of lenses with Ni mg > 
2, 3, 4, and 5 are quite similar. For example, the number 
of systems with Ni mg > 4 is about 1/3.5 of the number of 
systems with Ni mg > 2. 

It is interesting to ask for multiply-imaged quasars, 
what are the typical redshifts of lensing clusters and back- 
ground quasars. These are important for observationally 
identifying the lensing candidates. Fig.|6]shows the differen- 
tial number distribution for the lens redshift. The distribu- 
tions are similar in the ACDMO and WMAP3 cosmologies; 
both peak around redshift 0.5, and 90% of the lensing clus- 
ters are predicted to be between redshift 0.2 and 1. The pre- 
dicted lens red s hift di stribution is in good agreement with 
iHennawi etail <|2007bD . Notice that the median redshift for 
the lenses in the spectroscopic sample is slightly higher than 



Figure 7. The source redshift distribution of multiply-imaged 
quasars with separation > 10" in the ACDMO (left) and WMAP3 
(right) models. The peak of the source redshift distribution shifts 
to lower redshift in the photometric sample due to the lower mean 
quasar redshift for this sample (see f(3j) . The vertical bars show 
the sources redshifts for SDSS J1004+4112 (1.734) and SDSS 
J1029+2623 (2.197). 



that for the photometric sample. This can be understood as 
follows: quasars in the spectroscopic sample are on aver- 
age at higher redshift, since the most likely lens location is 
roughly mid-way between the lens and source, as a result 
the lens redshift is also slightly higher. The two observed 
cluster lenses have redshifts 0.55 and 0.68, roughly at the 
position of the peak, in good agreement with the theoretical 
expectations. 

The corresponding source redshift distribution is shown 
in Fig. [7] The distributions are again similar in the two cos- 
mologies. For the multiply-imaged quasars in the photomet- 
ric sample, their redshift distribution peaks around redshift 
2, while for the spectroscopic quasar sample, the peak is ex- 
tended with a range from 2 to 3. This is again a result of the 
lower source redshift in the photometric sample (see Sj3]) . The 
two observed source redshifts (1.734 for J1004+4112, and 
2.197 for J1029+2623) coincides with the peak predicted for 
the photometric sample, but appears to be at the low side for 
the the spectroscopic sample (but unlikely to be statistically 
significant due to the limited number statistics). 

It is interesting to explore whether these two cosmolo- 
gies are consistent with the data, given even the very lim- 
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ited number of cluster lenses. The effective area for the 
Data Release 5 of the SDSS is 5740 deg 2 . J1029+2623 
is from SDSS-II, which has now discovered about 10,000 
QSOs, which is roughly 1/4 of the QSOs discovered in the 
Data Release 3, which has an area of about 3700 deg 2 , 
so these two lenses were discovered in an effective area of 
~ 5740 + 3700/4 w 6500 deg 2 in the SDSS spectroscopic 
sample. In the ACDM0 model, the two crosses in the top 
left panel in Fig. [5] shows the positions we predict, (14.6, 
0.55) and (22.5, 0.93) for these two systems. The number of 
lenses with separation larger than 14.6" quadrupole-images 
(J1004+4112) and 22.5" double images (J1029+2623) are 
0.55 x 6500/5000 = 0.7 and 0.93 x 6500/5000 = 1.2. 
This is roughly consistent with what we have discovered 
in the SDSS so far. However, for the WMAP3 model, the 
positions of crosses in the bottom left panel in Fig. [5] 
are (14.6 0.072) and (22.5 0.094). The number of lenses 
with separation larger than 14.6" quadrupole-images and 
22.5" double images are 0.072 x 6500/5000 = 0.094 and 
0.11 x 6500/5000 = 0.14 for the SDSS spectroscopic sam- 
ple with an effective area of ~6500 deg 2 . We conserva- 
tively estimate the probability of observing such two sys- 
tems in WMAP3 model as follows. The predicted number 
of lensed quasar systems with Ni m& > 2 and 6 > 10" is 
0.36x6500/5000= 0.468. Assuming a Poisson distribution 
with a mean number of lenses A = 0.468, the probabil- 
ity of observing these two systems in WMAP3 model is, 
P(N > 2) = 1 - P(0) - P(l) = 1 - e- x - \e- x ~ 8%. So 
even with the very limited sample of two multiply-imaged 
quasars, it seems that the WMAP3 cosmology is compatible 
with observations only at ~ 8% level. Considering that our 
lensing optical depth may be over-estimated in the WMAP3 
model due to cosmic variance (see q4.4|l . and the fact that 
the number of observed systems may be a lower limit (due to 
possible incompleteness in the surveys), so the probability 
may be even lower (< 4%) in reality. 

4.3 Properties of central images 

For non-singular lensing potentials, catastrop he theory pre - 
dicts generically an odd number of images (|Burkd[l98lh . 
However, for nearly all g alaxy-scale lenses (with the excep- 
tion of PMN J1632-0033, IWinn et al.ll2003l . l2004h . the num- 
ber of images detected is even. It is thought that the central 
images may have been highly dem agnified by the stars at 
the centre of lensing galaxies (e.g., iNarasimha et all 1 19861 ; 
lKeetonll2003h . Furthermore the central images c an be swal- 
lowed by the presence of a c entral black hole ( e.g. , iMao et all 
l200ll ; iBowman et al.ll20o3 ; iRusin et al.ll2005l ). Even when a 
central image is present, it may not be easily identifiable 
because the central cluster galaxy is usually very bright. In 
this regard, colour information may be particularly useful in 
separating the (usually bluer) central QSO image from the 
central cluster galaxy. 

Extrapolating our experience for galaxy-scale lenses, we 
expect the central images on cluster-scale lenses will provide 
strong constraints on their central mass profiles. It is there- 
fore interesting to explore the properties of central images 
in cluster lenses. Below we discuss three image and five im- 
age cases in turn. Here, we take all of the images in the lens 
systems into account regardless how faint they are. For def- 
initeness, we choose a cluster population at z\ — 0.5 and a 



source population at z s — 2.5, roughly at the peak of the 
lens and source redshift distributions. 

The top panel of Fig. [5] shows the distribution of the 
differential cross-section of three-image systems in the plane 
of log /ii//i3 and R3, here {13 and /ti are the magnification of 
the faintest and brightest image, and R3 is the distance of 
the faintest image from the cluster centre. This plot shows 
two distinct, and disjoint regions. The bottom region arises 
due to three-image configurations illustrated in the top left 
panel of Fig. [3] in such cases, all three images are magnified 
and have similar intensity ratios. The region at much smaller 
R3 are due to three-image systems with a genuine central 
image close to the main cluster centre. 

In the bottom left panel of Fig. |8]we show the probabil- 
ity distribution of log/ii//i3. It can be seen most (~ 90%) 
of the three-image systems should have a central image 
brighter than 10% of the brightest images, and so should be 
relatively straightforward to detect. This conclusion, how- 
ever, ignores the two-image systems for which the central im- 
age has been infinitely demagnified. Such systems are shown 
as the hatched region, and is about 1/3 of the total number 
of three-image systems. If we take such systems as extreme 
cases of three-image systems with /ii//i3 — ► 00, the fraction 
will decrease to 0.9/(l+l/3)=65%. The probability distri- 
bution of R3 shows two prominent peaks (the bottom right 
panel of Fig. |HJ, one around 1" and the other around 10", 
followed by an extended tail. The peak around 1" is again 
due to genuine third, central images, while the peak around 
10" and the extended tail are due to systems as illustrated 
in the top left panel of Fig. [3] Our theoretical model appears 
to show for the two-image system, SDSS J1029+2623, the 
central image has a substantial chance of being observable 
given sufficiently deep images. 

Of particular interests are the central images in systems 
with at least four images, because for the first cluster quasar 
lens, J1004+4112, subsequent deep imaging reveals a faint 
central image, at abou t 0.3% of the brig htest image in the 
system (see Table 1 in llnada et al.|[2005l ). The top panel of 
Fig. [5] shows the number of five-image systems in the plane 
of log /ii //is , where /is is the magnification of the faintest 
image, and R$ is the distance of this image from the cluster 
centre. There are also five-image systems with R$ > 10"; 
they have been plotted in the right-edge pixels in the left 
panel. 

If we marginalise the separation from the cluster centre, 
we obtain the probability distribution of magnification ratio, 
/xi//i5, shown in the bottom left panel of Fig. [9] The solid 
line includes only those five-image systems with i?s < 5" 
while the dashed line includes all five-image systems. The 
solid line excludes most systems produced by merging sub- 
clusters. From these curves alone it seems that there is only 
a small probability of obtaining the magnification ratio of 
300, the observed value for J1004+4112. In fact, only 3% will 
have a central image be fainter than the brightest one by a 
factor of 100. However, this discounts those four-image sys- 
tems where the central image has disappeared. The hatched 
region is for those four image systems where the central im- 
age has been infinitely de-magnified by the BCG. The total 
cross-section for such systems is in fact quite large, about 
a factor of two of that for five-image systems. If we regard 
all these four-image systems as five-image systems with in- 
finitely demagnified central images, the fraction will increase 
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Figure 8. The top shows the normalised differential cross-section 
for three-image systems in the parameter space of log fii / fis and 
i?3 for the ACDMO cosmology, here R3 is the separation of the 
central image relative to the cluster centre, fii and /J3 are magni- 
fications for the brightest and faintest images; the latter is taken 
as the central image. The cluster is at z\ = 0.5 and the source is at 
z B = 2.5. The bottom left panel shows the differential probability 
distribution of log/^1/^3. The hatched region shows the corre- 
sponding cross-section for two-image systems, i.e., three-image 
systems with an infinitely de-magnified central image due to the 
BCG. The bottom right panel shows the probability distribution 
of R3. The cumulative distribution is shown as the dashed line. 

to ~ (0.03 + 2)/(l + 2) ~ 68%. If the profile is slightly shal- 
lower than the isothermal slope, p oc r~ 2 , these central im- 
ages will have a very small but finite magnification. Such 
central images may be what we are seeing in Jf 004+4112. 
Clearly the brightness of the central image will be a very 
sensitive probe of the central density profiles of BCGs. 

The bottom right panel of Fig. [5] shows the distribution 
of R5. It can be seen that, roughly 50 per cent of five-image 
systems have separation from the cluster centre smaller than 
1.5". About 25% have separation larger than 5", these are 
due to five- image systems produced by merging sub-clusters. 
If we include four-image systems, and assume all the in- 
finitely demagnified central images have R$ < 1.5", then 
the vast majority (~ (0.5 + 2)/(l + 2) ~ 80%) of these sys- 
tems will have central images within 1.5" of the BCG. These 
properties are broadly consistent with the observed system 
J1004+4112. 



4.4 Cosmic variance 

Our two simulations were run with a box size of 300/i _1 Mpc. 
A question naturally arises is whether they sample the 
mass functions well, especially for the high-mass end in 
the WMAP3 model, and how much our results are af- 
fected by the cosmic variance. The top panel in Fig. 1101 
shows the mass functions in our simulations. The solid 



Figure 9. The top panel shows the normalised differential cross- 
section for five-image systems in the parameter space of log fii/fis 
and R5 for the ACDMO cosmology, here R5 is the separation 
of the central image relative to the cluster centre, and fix and 
/Z5 are the magnifications of the brightest and faintest images 
respectively, the latter is taken as the central image. Systems 
with R5 > 10" have been collected at the right edge pixels. The 
bottom left panel shows the probability distribution of log fJ.i/p.5, 
which is simply obtained by integrating the left panel along the 
x-axis. The integration ranges are from to 5" for the dashed 
lines and from to 00 for the solid lines. The hatched region 
shows the corresponding cross-section for four-image systems, i.e., 
five-image systems with an infinitely de-magnified central image 
due to the BCG. The bottom right panel shows the probability 
distribution of R5. The cumulative distribution is shown as the 
dashed line. 



and dashed histograms are the mass functions at red- 
shift 0.5 in our ACDMO and WMAP3 simulations respec- 
tively. The solid and dashed curves are the theoretical 
ISheth fe Tormenl (2002) mass functions, which is based on 
the Press-Sch echter (1974) formalism and the ellipsoidal col- 
lapse model (|Sheth. Mo fc Tormenl l200ll ). We can see that 
the number density of halos in the WMAP3 model is always 
smaller than that in the ACDMO model on cluster mass 
scale. For M ~ 10 14 /i _1 Mq, the abundance of halos is lower 
by a factor of 3.1 in the WMAP3 model compared with that 
in the ACDMO model; for M ~ 10 15 /i _1 M©, the reduction 
factor is ~ 8.8. The number density in the last bin of the 
WMAP3 mass function is higher than prediction by a factor 
of ~ 6 due to the presence of just 6 clusters. The low num- 
ber of clusters at the high mass end in the WMAP3 model 
implies that ideally a larger simulation box is required to 
sample at least a similar number of massive clusters as in 
the ACDMO model. The two bottom panels show the dif- 
ferential cross-section as a function of the cluster mass in 
the two simulations. The total area under each curve is nor- 
malised to unity. The clusters are at redshift 0.5 and the 
sources are at redshift 2.5. 

The mass distribution of SDSS J1004+4112 was studied 
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Figure 10. Top panel: the solid and dashed histograms show 
the mass functions at redshift 0.5 in our ACDMO and WMAP3 
simulations respectively . The solid and dashed curves are the 
ISheth fc Tormenl ll2002fl mass functions. The two bottom panels 
show the differential cross-sections as a function of the logarithm 
of the cluster mass in the two models. The total area under each 
curve is normalised to unity. The clusters are at redshift 0.5 and 
the sources are at redshift 2.5. One can see that the box size of 
ACDMO simulation is large enough to sample the massive clusters 
but our WMAP3 simulation does not appear to have converged 
in the cross-section as a function of the cluster mass. 

in lOta et all (|2006h using X-ray data. Assuming isothermal- 
ity and an NFW profile for the cluster, the virial mass was 
derived to be M ~ 4.2±?;| x W^k^M®. The mass of SDSS 
J1029+2623 is more uncertain. The image separ ation implies 
a clu ster velocity dispersion a ~ 900km/s (see llnada et al.l 
120061 ) . and this can be converted to M ~ 5.6 x 10 14 /i _1 M© 
using cluster scaling relations, but with a large error bar. 
The box size of ACDMO simulation appears to be large 
enough to produce such massive clusters easily. On the 
other hand, there are only 20 clusters in the last three 
bins in the bottom-right panel of Fig. 1101 such a small 
number may be not sufficiently sample the clusters well 
in the mass-concentration space and result in a relatively 
large errors of differential cross-section at the very high-mass 
end. But the final optical depth is the integration of the 
mean cross-section as a function of redshift, which will sup- 
press the final variance. However our WMAP3 simulation 
does suffer seriously from the cosmic variance and appears 
to have too many very massive c l usters compared to the 
theoretical prediction. In iLi et "all l|2006bl ). we ran a lower 
resolution simulation in the WMAP3 model which evolves 
A?dm = 512 3 dark matter particles in a box with sidelength 
of 600ft -1 Mpc. Due to the larger volume, that simulation 
better samples the high mass tail of the cluster mass function 
and merger events. Comparing these simulations, we found 
that the smaller simulation over-estimates the number of gi- 
ant arcs by a factor of 2. For the same reason, the number 
of wide-separation lensed quasars in the WMAP3 may be 



over-estimated by a similar factor. The WMAP3 cosmology 
is compatible with observations only at ~ 4% and ~ 2% level 
for an over-estimation by a factor of 1.5 and 2 respectively. 
However, we caution that to more properly account for the 
cosmic variance, it is necessary to run a much larger simu- 
lation with high resolutions in order to sample the cluster 
mass function sufficiently. 



5 SUMMARY AND DISCUSSIONS 

In this paper, we have considered the number of multiply- 
imaged quasars lensed by clusters of galaxies in the SDSS 
photometric and spectroscop ic quasar samp l es. A s imilar, 
earlier attempt was made by IHennawi et al.1 l|2007bl ). This 
study extends their study by exploring how the predic- 
tions depend on the cosmology. We also study the multiply- 
imaged quasars as a function of image multiplicity, and ex- 
plored in some detail the properties of central images. Our 
main conclusions are as follows 

(i) We found that the predicted multiply-imaged quasars 
with separation > 10" in the SDSS photometric sample 
(with an effective area 8000 deg 2 ) is about 6.2, and about 
2.6 in the spectroscopic sample (with an effective area 5000 
deg 2 ) in ACDMO model. These numbers are reduced by a 
factor of ~7 or more in the WMAP3 model. 

(ii) The predicted cluster lens peaks around redshift 0.5, 
and 90% are between 0.2 and 1. This distribution is largely 
independent of cosmology, and similar for both the photo- 
metric and spectroscopic samples. 

(iii) The relative number of systems with Ni mg > 4 im- 
ages and those with Ni mg > 2 images is about 1/3.5. This 
ratio is largely cosmology independent, and similar for both 
the photometric and spectroscopic quasar samples. 

(iv) Because we modelled the BCGs as a truncated 
isothermal sphere, this creates a region, comparable in size 
to the angular Einstein radius of the BCG, inside which the 
central image disappears. For most central images in five- 
image configurations are quite faint, and close to the clus- 
ter centres. For three-image systems, the central images are 
brighter and further away from the cluster centres. 

One uncertainty in our calculation is the ad hoc inclu- 
sion of a BCG at the centres of clusters at all redshifts. In 
reality, the BCGs are assembled as a function of time, and so 
it is interesting to estimate the effects of an evolving popu- 
lation of BCGs. We adopt a simple model by assuming there 
are no BCGs at the centres of clusters above redshift 0.5. 
Hennawi et al. (2007) finds that the optical depth decreases 
by a factor of 1/3 if all the BCGs are ignored. However, from 
Fig. 6, the contribution of clusters with Z\ > 0.5 is about 1/2 
of the total cross-section. Thus if we remove all the BCGs 
from the clusters above redshift 0.5, the optical depth will 
decrease roughly by ~ 16%. So the evolution effect may be 
modest. In any case, such a decrease in the optical depth 
will make the match between the WMAP3 model and the 
observations slightly worse. 

Our predicted numbers of multiply-imaged quasars for 
both the pho tometric and spectrosc opic samples are lower 
than those bv lHennawi et all l|2007bl ) roughly by a factor of 
2 in the ACDMO cosmology. Part of this diffe rence is due 
to the slightly higher ag (0.95) adopted by IHennawi et al.l 
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(2007b). The abundance for typical lensing clusters around 
M > 10 14 - 5 /i _1 M© at redshift 0.5 is about 30% higher than 
that in our model with as = 0.9. Their higher as not only 
provides more massive clusters but also makes the clus- 
ter formation earlier and with higher concentration; both 
may increase the predicted number of cluster lenses. Our 
quasar lumino s ity fun ction is similar to those adopted in 
iHennawi et al.l (|2007bl ) for the low redshift intervals, but is 
a factor of 2 higher than that adopted in their paper for 
the high- redshift interval (z 3 ~ 4). So the difference in the 
source luminosity function goes in the opposite direction. 
Note, however, the number of such higher redshift quasars 
is small and so this should not impact significantly on the 
overall predicted large-separation lenses (see Fig. [TJl . 

For galaxy-scale lenses, the number of lenses with 
Nimg > 4 is about one half of that with Mmg > 2. For 
example, the statistical sample of the Cosmic Lens All Sky 
Survey (CLASS) has 13 le nses, 6 have at least 4 images (in - 



cluding one with 6 images, iBrowne et al"1l2003l ; ICIiaell2003f ). 
Our study indicates the ratio is about 0.3 for cluster-scale 
lenses, about a f a ctor o f 1.5 smaller. This is consistent with 
lOguri fc Keetonl l|2004h . Using a triaxial cluster model, they 
found the ratio is 0.2 for inner density profile a = 1 and 
0.4 for a — 1.5. The comparison between the two studies is, 
however, not straightforward as we use different definitions 
of iVi m g and their study also all of the faintest images in 
the multiply imaged systems can be detected, while in our 
study we assume certain magnitude limits for the images 
to be observable (see H3.1|l . If we relax the magnitude limit 
to i — 25, then we find that our ratio is changed to 0.35, 
still w ithin the bracket of the two values of lOguri fc Keetonl 
il2004h for differen t inne r den sity slopes. As discuss ed by 
iRusin fc Tegmark I (|200ll ) and lQguri fc Keetonl (|2004h . this 
ratio depends on a number of parameters, including the 
shape of the lensing potential and the density slopes. It ap- 
pears that this is mainly because of the shallower inner den- 
sity slope which causes the factor of 1.5 difference in the 
relative numbers between cluster lenses and galaxy lenses. 

One striking result of this study is the reduction in the 
predicted number of large-separation lenses in the WMAP3 
model compared with that in the ACDM0 model. This re- 
duction is comparable to that for the giant arcs (~ 6) in 
these two cosmologies. However, we emphasise that the un- 
certainty for the quasar sample is smaller, because the source 
population is reasonably well known (see (J3J. There are two 
lens systems, J0114+4112 and J1029+2623, have been dis- 
covered in the SDSS spectroscopic sample with an effec- 
tive area of ~ 6500deg 2 . In the ACDM0 model, from Fig. 
[S] the number of lenses with separation larger than 14.6" 
quadrupole-images and 22.5" double images are 0.7 and 1.2. 
This is roughly consistent with what we have discovered in 
the SDSS. However, for the WMAP3 model, the number of 
lenses with separation larger than 14.6" quadrupole-images 
and 22.5" double images are 0.09 and 0.14 for the whole 
spectroscopic sample. Assuming a Poisson distribution, the 
probability to observe two quasar systems in the WMAP3 
model is ~ 8% and decreases to 4% or lower if we account 
for the cosmic variance in our simulation (see H4.4p . So it 
appears that a higher as or a higher Qm,o is preferred by 
the data. 

Clearly it is desirable to have a much bigger sample of 
large-separation multiply-imaged quasars. Our study shows 



that most cluster lenses have already been discovered in the 
SDSS spectroscopic quasar sample, but several more candi- 
dates may yet to be discovered in the photometric sample. 
Future large surveys m ay discover many more such examples 
|Wittman et al.ll2006l ). which will provide strong constraints 
on the cosmology and the central mass profiles of clusters of 
galaxies. 
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